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Abstract 


Many real-world problems occur under uncertainty. In this paper, we 
consider interval linear programming (ILP) which can be used to tackle un- 
certainties. Several methods have been proposed by researchers, such as 
the best and worst cases, Two-step method (TSM), improved TSM, ILP, im- 
proved ILP, three-step method, and robust two-step method. First, we define 
feasibility and optimality conditions in ILP models and review some solving 
methods shortly, and then show that some solutions of the TSM method 
are not feasible. Therefore, we propose an updated TSM method (namely, 
UTSM) by considering the feasibility and optimality conditions. In this pa- 
per, the UTSM method was applied to identify the reduction of aerosols by 
using two controllers with a minimized cost to demonstrate its application 
under uncertainty. Compared with other methods, the solutions obtained 
through ILP were presented as interval, which can provide intervals for the 
decision variables, objective function, and decision-makers. Therefore, the 
decision-makers can make the best decision based on the obtained solutions 
through ILP, and then identify desired plans for aerosol-emission control un- 
der uncertainty. 
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1 Introduction 


Interval linear programming (ILP) is described as the lower and upper 
bounds, which is used as a powerful tool to deal with optimization prob- 
lems under uncertainty. Some researchers have been working on solving 
methods of ILP. In the best and worst cases (BWC) method proposed by 
Tong [80], the ILP model was converted into two submodels: the best and 
worst submodels, which they, respectively, have the largest and the smallest 
feasible regions. A given point is feasible for the ILP model if it satisfies the 
constraints of the best problem, and it is optimal for the ILP model if it is op- 
timal for at least one characteristic model. The BWC method was extended 
by Chinneck and Ramadan for ILP models with equality constraints; see [7]. 
A novel ILP method was proposed by Huang and Moore [18], and Huang and 
Cao [17] also proposed their analyzed principals of two-step method (TSM). 
Some solutions obtained through the BWC, ILP, and TSM may be infeasible. 
New solving methods named three-step method (ThSM) and robust two-step 
method (RTSM) have been developed for solving ILP models; see [10, 17]. 
To guarantee that the given solutions of ILP method are completely feasible, 
Zhou et al. [33] proposed the modified ILP method (MILP) in which an ex- 
tra constraint is added to the second submodel. Since the resulting solutions 
from the MILP may be nonoptimal, two improved ILP and MILP (IILP and 
IMILP) methods for solving ILP problems have been proposed; see [4]. The 
solutions of these methods are completely feasible and optimal; see [3]. Also, 
there are other methods for solving ILP models including IThSM, a new al- 
gorithm by Ashayerinasab et al. and a method proposed by Garajov et al.; 
see [2,5, 12). 

In this paper, we update the TSM (namely, UTSM) by considering the 
feasibility and optimality conditions. Following that, a numerical example is 
solved to demonstrate the effectiveness of UTSM. The air pollution caused 
by polluted resources, which cause pollution emissions in many sources, re- 
ceptors, and artificial factors, including emissions of dangerous gases from 
factories and coal-fired power plants which pose serious threats to humans. 
These can have a greater impact under stable atmospheric conditions and in- 
creased pollution emissions. Pollution control and the reduction in the levels 
of atmospheric pollution are very important. These would approach the lev- 
els of atmospheric stability and prevent the destructive effects of pollutants. 
The proposed models can identify pollution sources and atmospheric weather 
conditions. These methods suggest pollution control using efficient methods 
for the study area. In this case, the levels of emissions of pollutants and 
environmental loading capacity need to be considered. Using optimization 
models, the costs should be reduced by controlling methods. 

Another aim of the study is to check air pollution management under 
uncertainty as a range in the ILP model. The suggested methods will be 
analyzed to show the control and the reduction of dust by mulching and 
green belt controllers. However, the suggested methods are common methods 
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for dust and aerosol control. Now, to assess the effectiveness of the two 
controllers (mulching and green belts), the ILP model will be used and solved 
by the BWC, ITSM, and UTSM methods. The aim of this model is to 
minimize the costs in order to reduce aerosols by using two controllers. 


2 Interval linear programming 


In this section, we define the ILP model and some solving methods. An 
interval number [X~, Xt] is shown as X* where X~ < Xt. If X7~ = X7, 
then X* is degenerate. If A~ and A* are two matrices in R™*” such that 
A~ < At™, then the set of matrices At = [A~, At] = {A|A~ < A< AT} 
is called an interval matrix, and the matrices A~ and A® are called bounds 


of A. Center and radius matrices are defined as Ays = 5 (A*—A7) and 
A® = 3(A~ + At). A special case of an interval matrix is an interval vector 
x* = {x|x~ <x < x7}, where x ,xt € R"; see [1]. If A* is a square interval 


matrix and each A € A* is nonsingular, then A* is called regular. Consider 
the following ILP model: 


n 

S.t. > AjjX5 < bi, = LD egy 
j=l 

x; > 0, ee ey ere LF 


where ajj € a 
is defined as 


ij¢7 € C7, and b; € b;. The feasible solution set of the ILP 


n 
ze R”: Se Se leat a le 
j=l 


Also, the optimal solution set of the ILP is defined as the set of all optimal 
solutions over all the characteristic models. Some methods for solving ILP 
models are BWC, TSM, ITSM methods (see [16,30,31]), which are presented 
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in Table 1. Also f;,, and an are the best and worst optimal values of the 


objective function and the optimal solution is 27,4 = [aj opt» Tp opt]: 


We also have 


la] 7 ai; Z 0, | + )+ an ee ai; 2 0, 


—a;,, aj; <0, 


Note that in all methods, for the objective functions, the notations “++” 
and “-” are always used for the first submodel and the second submodel, 
respectively. But for the variables, in the TSM and ITSM methods, the 
notation “+” for 7;, 7 =1,..., k, in the first submodel and 7 =k+1,...,n 
in the second submodel is used. Also, the notation “-” for 2}, 7 =k+1,...,n 
in the first submodel and j = 1,..., k in the second submodel is used. In 
the BWC method, these notations are not used. 


Remark: If a decision-maker wants to minimize the costs, then the BWC 
method can provide the lower and upper bounds of the cost. Now, an exam- 
ple is solved and then the results are compared. Consider the ILP model as 
follows: 


8.t. 
[1.1, 1.2] ¢F + [1.7, 1.9] a> < [11.7, 12.1], (2) 
[4, 5] ty [3, 4 ra < [5, 7 ’ 
O<a, <a], 
O<2% < 25 


The results obtained through the BWC, TSM, and ITSM are given in 
Tables 2 and 3 and Figure 1. 


Some points of the solution space obtained by the BWC and TSM are 
infeasible. For example, the point (5.3839, 4.0076) in the BWC and the 
point (5.1417, 4.3388) in the TSM are infeasible because they do not satisfy 
the inequality 1.laf + 1.7ry < 12.1, which is the first constraint of the best 
problem. The solution space obtained by the ITSM is feasible. 


In the next subsection, we will define optimality condition and show that 
some points of the BWC, TSM, and ITSM are nonoptimal, and then we 
propose a new method, which we call an updated TSM method. 
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submodels of BWC, TSM, and ITSM 


Table 1: 


u ‘ oe ane a 
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ui ‘tT +y=l' (< '5)uBrg e -e)UsIg {2 a 
eee ae od a Sayaaeg = (Seyutig J Peat uc Tty=l og f 
et fo)u81g (ze)uais f2 ae (in yo TSel G22 
=f Sta — (! 3 ydo £ t+y=0 T=o 
ps fae | ye '5)u8tg (2 Fe)usig f2 Ee la as Le Seas [te <q rae Fx( Neyusig _ [fe] ~< 
ot—C alto d 4 WS.LI 
49> 5 82 <4 ae 
1 ty cc eer ts! Fn ct a Sd 
‘qs ae 2e)UsIS | fre] << + =x(¢e)usIg , | fn| < Ja fo << + jejo < =f xeu 
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Table 2: submodels of BWC, TSM, and ITSM for ILP model (2). 
Methods The first submodel The second submodel 
max ft = 427, — 1.529 max f = 3.5%; —1.7x2 
s.t. S.t. 
4a, — 429 <7 5x21 — 3842 <5 
%1,%2>0 %1,%2>0 
max f  =3.52%, —1.7a3 
max ft =4a} — 1.523 s.t. 
get. 12a; + 1.72f < 11.7 
TSM Lia} +1927 < 12.1 5a, —3af <5 
Agi —4a5 <7 #y < aio 72] < 5.1417 
efi, =O €3 > Cy ope > 13 > 3.3917 
Xi oe >0 
max f— =3.5a,7 —1.7a7 
S.t. 
max ft = 4a} — 1.523 1.227 +1.7af < 11.7 
8.t. 5a, —3at <5 
ITSM lla} +1927 < 12.1 Lda pet Lteg << 12.1 
4gt —4ry <7 Att oot — 402 <7 
a], ay > 0 ty Seo 
iis 2 L2 opt 
ta 0 


Table 3: Solutions of BWC, TSM, ITSM for ILP model (2). 


Methods x7 La 2* 
BWC [3.4046,5.3839] | [3.6339,4.0076] | [5.1031,16.0848] 
TSM (3.6033,5.1417] | [3.3917,4.3388] | [5.2355,15.4792] 


ITSM 


[3.2744,5.1417| 


[3.3917,3.7906] 


[5.0162,15.4792} 
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—— The best and the worst constraints "=" TSM 


—. BWC om ITSM 
6 


xl 


Figure 1: Solutions of ILP model (2) 


2.1 A new method: An updated TSM 


In this subsection, we suggest an updated TSM by considering both feasibility 
and optimality conditions. At first we review basis stability, and then obtain 
two feasibility and optimality conditions. 


Definition 1 (see [15]). The problem max {c’x: Ax =b, x > 0}, wherec€ 
C+ CR”, Ac A* CR™*” and b € b= CR” is called B-stable with basis B, 
if B is an optimal basis for each characteristic model. The ILP model is called 
(unique) nondegenerate B-stable if each characteristic model has a (unique) 
nondegenerate optimal basic solution with the basis B. Let B C {1,2,..., n} 
be an index set such that (the restriction of A to the columns indexed by 
B) is nonsingular. Similarly, N = {1,2,...,n}\B denotes the index set of 
nonbasic variables and Ay is the restriction of A to the columns indexed by 
N. Also, B can be computed by solving an arbitrary characteristic model. 
We now review the conditions for B-stability. 

Regularity: Ag is regular. 

Feasibility: The solutions set of the interval system Apxp = b are nonneg- 
ative. 

Optimality: Ap is optimal, that is, chA, An —ch >07. 


Theorem 1. [27] If p ((45" Aas) <1, then Ag is regular, where p(.) 
denotes the spectral radius. 
Theorem 2. /27/] If maxi<i<n ((405"| Aas) > 1, then Ap is not regu- 


lar. 
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Some conditions for the regularity of interval matrices have been proposed 
in [29]. 


Theorem 3. [28] The interval vector r is an enclosure to the solution set 
of a system Apup = b, where 

= min {—2f + (ef + |efl) Mis) + aye (—2t + (36 + [ofl Ma} 

rf = max {af + (@f — |ael) Mis) + sab (et + (28 - los) Mu, 
M=(I- a Raga = (AS), 

xt = M(\2*| + |(4°)5'| An), 


r, 


with AG is nonsingular and p ((405"| Aas) <i. 


Theorem 4. /15/] Let y be an enclosure to the solution set of Apy = cg. If 
((AR) y) > Gi, then the optimality condition holds. 


Theorem 5. /15/] Let diag(q) denote the diagonal matrix with entries 
Q1;+-+;Qm- For each q € {+1}", if the solution set of the system 


((A%)* — (Aa,)” diag(q))y < 
—((Ag)* + (Aan) diag(q)) y 
diag(q)y = 9, 


cs 
< i 


lies in the solution set of the system 


{a —(Aay)’ diag(a)) y < ch, 
diag(q)y = 0, 


then the optimality condition holds. 


Theorem 6. /15/] Let the ILP model be unique, nondegenerate B-stable, 
where B = (1,...,m). Then the optimal solution set of the ILP model is 
equal to the solution set of the interval system 


Aptp =b, xp >0, tn = 0, or At rp 2b, Al@B < b+, xp >0, ey = 0. 


Now, we introduce two submodels of UTSM. Firstly, we solve the sub- 
model corresponding to f~ ,which is the second submodel of the TSM. Sec- 
ondly, we add two extra constraints in order to ensure that the resulting solu- 
tion space is completely feasible and optimal in the submodel corresponding 
to f+. Submodel 1 is as follows. 

Submodel 1 
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max f— -39 a, + 3 Gay 


j=k+1 
8.t. 
k 1 
> Jaxj|* Sign (az; )xyp + 3 lai] Sign(a;)x;" <b>, i=1,... 
j=l j=k+1 
z, 20, j=1,...,h, 
a; >0, j=k+1, >in, 


For Submodel 2, x5, for j =1,..., k and z 
optimal solutions of the first sibmiadel. 


Cota = Pda, 


j=l i j=k+1 
8.0. 
k 
> laig|” Sign(ai x + 3 lay|* Sign(az)x> < bt, i=1,. 
j=l j=k4+1 
OS 2 opt aa PS Nady hy 
O< 2; aes j=k4+l,...,n, 


Now, we obtain an extra constraint to ensure that the solutions 
mal. 


Theorem 6 implies that » Oe jtj 2; , t= 1... 


k n 
ao ae = 
p31 ata; 7 > ay; C i> ps] As L i : > apa; ;2 b; 
—pt+l j=k+1 j=n—qtl 


Since ai; See =1,...,k-p, 


a 


ax <0 forj=k—-ptl.,..., 
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(3) 


n are the 


are opti- 


1,...,n—qgand aj, >0 forj=n—q-+1,...,n, therefore 
k—p — 
S a, ij + » a; Lit » ayj;Cj> » Qi ij 
I= j=k—ptl Gak+ oe 
= ae “ + xt ae 
2 a, ih; opt b » Ajit a = oi er 3 a, ie; . 
j=l j=k—-p+l1 j=k j=n—qtl1 
Therefore, for ee it is sufficient iat 
k—p 
aja, gj opt + 3 a, am > a OF ee >a ai, Jj Zz b; ’ 
j=l j=k—p+1 j=n—qtl 
or 
vio vt tf Sign(a;) = Sign(cy), j=1,...,k, 
a Lf Sign(as Sign(c*), j=1,...,k, 
S ae! ss be v= ei if ign( ij) Fa ign( i) J 
ja Tron tf Sign(aj) = Sign(c"), j=k+1,...,n 
x; if Sign(az) # Sign(ct), j=k+1,. 


To ensure that the solutions are feasible, let us obtain an extra constraint: 
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Theorem 6 implies that s a,,05 < Bh Nada gO 
j=l 
k—p k n—-q n 
ati + DY ayajs+ DY agajt Dd) ajay < OS Lyn ss 
j= j=k—p+1 j=k+1 j=n—qtl 
Since ai; > 0 forj =1,..., k—p, a; <0 for j = k—-p+1 K+1,...,n-@ 
and 
aj; 2 0 for j =n—-—q+1,..., n, therefore 
k—p n—q n 
yee + DP ageit DD ayayt DL ay48; 
j=l j=k-ptl j=k+1 j=n—qtl1 
k—p exc k n-q a 
< xX aja} + » ai; Tjopet x ajja; + > i525 ont 
=1 j=k—-pt+l1 j=n—qtl 
Therstate. for feasibility, it is sufliciént: that 
k—p k n—-q n 
y= ayia} + > Qi 5L 5 opt t ajjt; + > QiZLj opt S by or 
j=l j=k—pt j=k+l1 j=n—qtl1 
ay if Sign(a.) =Sigen(cs),. 7= 1,..<,1, 
3 a,2’. wee br, ql = 5 opt if Sign(ajj) Sign(cy), . j —_ 1, . ail 
= i Zl if Sign(aj) = Sign(c;"), j=1+1,...,0, 
Li opt uf Sign(a;;) # Sign(c"), j=1+1,...,n, 
Finally, we have fon, = [ Fowtyr faves West = (Beers pee] + “Theretore, 
the second submodel is as follows. 
Submodel 2 
+_oa tet - 
max ft = Di cpap + Di fa; 
j=l j=k 
st. Y> |ayj| Sign(aF)xF + = |aij|* Sign (a3 ay )x; < bt, i=1,...m, 
j=l j=k+1 
OS 2 opt < ay, a eee 
O52; Soop» J=R+I, i, — 
© opt ene a ee ne j=l, _ 
nm ee aye Sp Se Lines 
ate" > b-, gl. = a v ign(a ij / IgNn{c ji? J ’ >™) 
p> = : Xs opt if Sign(a ai) = — Sign(c cj), J — k+ 1 
x, tf Sign(a) A Sign(c*), j=k+1,... 
Tope tf Sign(az) = Sign(ce), j=1,.-..1 
Oo x; if Sign(a;) A Sign(c*), j=1,...,1 
— a! < bt ae J ij / ji? ’ at) 
2, ou? SOP, Trop: tf Sign(aj) =Sign(c), j=l+1,...,n 
a if Sign(a;;) 4 Sign(c;"), j=1+1,n. 


+ yl! 
jjej; 2 6; are feasibil- 


n 
Therefore, the constraints x aya,’ < bf and x a; 
=1 


ity and optimality conditions, creapactivalie 


Note that, since the solution regions obtained through the UTSM method 


is completely optimal, the union of these regions will be closer to the exact 
optimal solution region of the ILP model. 
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Now, consider ILP model (2), again. Given Figure 1, the points 
(3.4046, 3.6330), (3.6033, 3.3917), and (3.2744, 3.3917) obtained by the BWC, 
TSM, and ITSM respectively, are nonoptimal, because they do not satisfy 
inequality 1.27;+1.9x2 > 11.7, which is the first constraint of the worst prob- 
lem with the reverse sign. Model (2) is equivalent to the following model: 


max z+ = (3.5, 4Jaf — [1.5, 1.7] 25> 


8.0. 
[td L2)ar + (7, 18) es +as = (11.7, 12.1), 
[4, 5] ry _ [3, 4] ao) . v4 = [5, 7) 
<2, <<, 
<a, <a, 
<a, < 27, 


Oia, A a7. 


According to Theorems 1, 3, and 4, ILP model 3 is B-stable, where B = 
(1,2). Now, we solve the model by the UTSM. 
submodel 1 
max z~ =3.52, — 1.72} 
s.t. 
1.Q¢7 +1.?72f < 11.7, 
5a; — 3a¢ <5, 
“, te > 0, 
submodel 2 
max z+ = 4a} — 1.527 
s.t. 
Lie} +199; < 12.1, 
4ay —4Arg <7, 
Ly > Ly op = 3.6033, 
Ly < oi = 4.3388, 
Lley + Le dees S 12.1, 
ast —4ry <7, 
1 Dig 1.1 ons 
{ bar — 3x, 2 5X1 opt = SUF oes 
a7, 


Now, the results obtained by the UTSM have been given in Table 4 and 
Figure 2. 


Table 4: Solutions of UTSM for ILP model (2). 


Methods Lp xy = 
UTSM [3.6033,4.2945] (3.8820, 4.3388] (5.2355,11.3550] 


As shown in Figure 2, the UTSM introduces a feasible and optimal space. All 
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points satisfy the constraints 1.la, + 1.7% < 12.1 and 4a, — 4a < 7 (feasi- 
bility conditions) and the constraints 1.27%; +1.9%2 > 11.7 and 5x, — 3a > 5 
(optimality conditions). 


—— The best and the worst constraints UTSM 


6 


Figure 2: The solution region obtained by the UTSM for ILP model (2) 


3 Application to air quality management 


Atmospheric conditions can be affected by many factors, including SO2, CO, 
aerosols, dust, and so on. Aerosols have a very high diversity. As one of the 
greatest pollutants, this phenomenon occurs in arid, semiarid regions, and 
deserts, such that their sizes are different according to the concentrations in 
air. These can have destructive effects on health, nature, and the environ- 
ment and create numerous economic problems for local people. It can be said 
that sandstorms and dust-storms are natural events with significant concen- 
trations of particulate matter, which depend on the wind speeds. When the 
wind speed is almost more than 8 meters per second (and sometimes based 
on the region, which is stable or not, the wind speed is more than 6 me- 
ters/seconds), it can get into the air flow and the earth’s atmosphere based 
on the roughness, particulate size, soil texture, soil moisture, and vegetation; 
see [9,32]. The suggested methods are mulching and green belt that should 
be analyzed to control and reduce dust. These suggested methods are con- 
ventional methods for dust control, and many studies have been conducted 
in this area. 


An updated two-step method for solving interval linear programming ... 219 


Mulch is a non-life coverage used for sand dune fixation in the desert 
lands. A variety of organic and inorganic mulches are available that can be 
used to accelerate production, reduce erosion, and even help control weeds. 
In Iran, oil mulch is usually used to cover. Most of the causes of oil mulching 
are the immediate effect, using them is provided in a short time, material 
source is available in the country, and there is no need for foreign resources. 

Vegetation can have a major impact on air pollution and climate quality. 
It can also directly or indirectly decrease the temperature, increase moisture 
and dust absorption, and affect urban air pollution control. Studies on this 
subject have concluded that trees can greatly control air pollution, prevent 
harmful pollutants like PM10, SO2, NO2 and CO, and particulate matter 
in the atmosphere. Studies have shown that a lot of PM10 and PM2.5 is 
destroyed by trees, such that the efficiency and effectiveness of this factor 
plays a very important role in social and human health, and air quality 
control. In many countries nowadays, trees are being planted to reduce costs 
and harmful factors in the atmosphere; see [19, 21, 22]. 


Figure 3: Trees as air quality buffer. 


Particulate 
pollution 


— 


Figure 4: Trees as air quality buffer. 


It has been well-documented that the plants destroy air pollutants, in- 
cluding hydrogen fluoride, sulfur dioxide, and some combinations of photo- 
chemical reactions in the air. It can be said that the green belt efficiency in 
the reduction of pollution defined by the following formula: 
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As F is the relative reduction effect by trees (percent), R is amount re- 
moved by trees (kg), and A is the amount of pollution in the atmosphere (kg). 


3.1 Methodology and overview of the study system 


The Sistan and Baluchestan region in Iran is located in the southeast, with 
arid and semi-arid weather. The Sistan plain covers a wide range of the 
western part of the province, and almost a wide catchment area. In addition 
to Iran, it has expanded into two countries Afghanistan and Pakistan. The 
border areas between Iran, Pakistan, and Afghanistan are the main dust 
source region in southwest of Asia; generating some 81 dust storms over 
Sistan; see [13, 24-26]. The Sistan region is one of the richest regions of 
Tran in terms of wind, and this sometimes creates difficulties for its people. 
Dust storms generate great waves, especially in Sistan. Figure 6 shows the 
geographical location of the Sistan and Baluchestan region. 


1. 


Pakistan 


Balochistan | 


ie 


€2010 Google - imagery ©2010 TerraMtetrics, Map data ©2010 LeadDag Consulting. AND. ? 


Figure 5: Geographical location of the Sistan and Baluchestan region 


The Sistan’s 120-day winds are the most famous local winds of Iran, 
blowing from mid-May to mid-September, over a large area in the northern 
part of Sistan-Baluchestan. The strong “Levar” winds in summer favor the 
uplift of large quantities of dust from the Hamoun basin, which is located in 
the northern part of Sistan; see [25]. The 120-day winds of Sistan-in terms 
of time and local distribution and wind speed-have a high variety. In Sistan, 
the strong winds above eight nodes start in May and continue till November. 
In Sistan, winds are severe in summer, but can be slowed down in other 
seasons. Therefore, in this section, we focus on the optimization of an air 
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quality management model and PM10 emissions control in the Sistan region 
in a source through controllers, as mentioned in the previous section. 

The aim of an air quality management model is to control effects on health 
by air pollution. These effects are directly related to pollutant concentration. 
Now, we establish the relationship between air pollution emissions and load- 
ing concentrations. The ground concentration at each downwind location 
(x,y) can be estimated as follows [14]: 


Q <2 H?2 
ee TUT yO z maid (=) i (s3)| , 


where C'(a,y) is the ground-level concentration of pollutants (mg/m3) at 
point (a, y); x is the downwind distance from the source; Q denotes pollu- 
tant emission rate (mg/sec); u is average wind speed (m/ sec) along x direc- 
tion; and oy and a, represent horizontal and vertical dispersion coefficients of 
plume versus(against) downwind distance x from the pollution source. Gen- 
erally, 0, and a, can be approximated by the following equations [8, 23]: 


Oy = ya >, Ge= yor? 


Wind Direction /. ub (Up 


Figure 6: Dust-plume geometry 


where 7; and /; are constants along the y direction and yz and {2 are 
the coefficients along the z direction, respectively. The values of 71, 81, ye, 
and (2 are related to atmospheric stability classes defined for different me- 
teorological situations through wind speed, solar radiation (during the day), 
and cloud cover during the night [14, 23]. 
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Figure 7: Typical visibility within a dust-storm 


A transfer factor t, can then be determined as follows [14]: 


P TUGyOx ‘ 202 Qo2) |’ 


where t, represents the contributions of pollutant emission rate at the Sistan 
area to the ground concentration of receptor zone p; the index p indicates 
different receptors influenced by air pollution emissions; H indicates the ef- 
fective height of the air pollution plume from the Sistan region [6]. 


3.2 ILP model for air quality management modeling 


The optimization model of air quality management has been presented as 
fuzzy linear programming [11,20]. Now, we consider the following ILP model: 
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where 7 is the PM10 emission source, j is the type of PM10 control method 
(j = 1, 2, respectively, denotes mulching and green belt), and & is the plan- 
ning period. The planning horizon is five years. Decision variables are the 
amount of PM10 allocated to control measure j for the reduction from source 
i in the period k (t/day). It has been shown by the symbol X;;,. The ob- 
jective is to minimize the total cost for the abatement of PM10. Now, the 
model describes the study area according to particulate reduction methods 
in this article. 
Also L;, is length of period k (day), T'Cj, is operating cost of control mea- 
sure j during period k ( Rials/t), FCj, is maximum mitigation capacity of 
measure 7 in period k, ? is PM10 generation amount of emission source in pe- 
riod k (tonne/day); 7; is efficiency of control measure j, ex is PM10-emission 
allowance for Sistan during period k (tonne/day), t, is transfer factor from 
emission area to receptor zone p (day/m3), and 6;, is environmental loading 
capacity of receptor zone during period (mg/m3). Note 7 = 1 represents the 
emission source of PM10 in the Sistan region. All the parameters in model 
(4) are presented as interval numbers. In this case, the parameters of TCjx, 
FC, Sk, j, €r, and Oxy are assumed to be estimated as interval numbers. 
Particulate matter in the Sistan region is discussed over four-month periods 
(the maximum amount of dust and particulate matter) during five years. We 
estimate-according to two methods, mulching and green belts-the tables are 
related to the costs of each method, efficiencies, and permitted amounts of 
dust in the region, and show their reduced values by using the two methods in 
the respective tables. Here, S; is the aerosol generation amount of emission 
sources in period k in Sistan, which has been shown in Table 5. 


The parameters related to the model in the Sistan region are now de- 
scribed. T'c;, represents mulching costs per hectare with respect to the in- 
flation rate (million). If we require 10 tonnes mulchs for every hectare and 
is equal to 20,000 Rials per liter of mulch, then the amount calculated using 
the inflation rate will be equal to Table 6. 
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Table 5: S;-emission rate for ILP model (4) 

Months 2012 2013 2014 2015 2016 
May [55,1345] (76.6,4985] | [62.5,4983.6] | [15.3,460.1] | [45.6,4985] 
June | [130.7,1384.8] | [60.1,2570.4) | [140.2,4985] | [17.2,452.4) | [49.1,630.8] 
July [315.6,9985] | [64.6,3044.6] | [4979.4,4985] | [10.5,4985] | [56.2,2489.4] 

August | [100,1455.7] | [280.2,3553.5] | (130.2, 4985] | [880.9,4840] | [92.2,4983.6] 

Table 6: Costs of different pollution controls for ILP model (4) 
T Cir 2012 2013 2014 2015 2016 
Mulching | [8.9,18] | [12.1,24.1] | [15,30] | [16.8,33.8] | [16.6,37.3] 
Green belt | [2.5,6] 3.6] B,7] 4,3] [5,10] 


Figure 8: Dust precipitation 


where e;, denotes the aerosol (PM10) emission allowance for the source 


(Sistan) during period k (t/day) . 


If FCjx is the maximum mitigation ca- 


pacity of measure 7 at area in period k, then the amount of the reduction by 
the mulching and green belts methods can be seen in Tables 8 and 9. 


Table 7: e, for ILP model (4) 


PM10-emission allowance PM10 
Good 0-54 
Average (allowance) 55-154 
Unhealthy for sensitive groups 155-254 
Unhealthy 255-354 
Very Unhealthy (crisis) 355-424 
Dangerous (emergency) 425-604 
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Table 8: The value of PM10 emissions controlled by mulching 


Months 2012 2013 2014 2015 2016 
May [44,1076] [61.3,3988] [50,3986.9] [12.24,368.08] [36.48,3988] 
June [L04. 6,1108] [48.08,2056.3] [112.2,3988] [13.76,361.9] [39.28,504.6] 
July [252.5,7988] [51.68,2435.9] [4979.4,3989.1] [8.4,3988] [44.96, 1991.5] 

August [80,1164.6] [224.16, 2842.8] (104.16, 3988] [704.72,3872] [73.76,3986.9] 
FCjk [120,2834.1] [96.3,2830.7] [1311.4,3988] [184.8,2147.5] [48.62,2617.8] 


Table 9: The value of PM10 emissions controlled by green belt 


Months 2012 2013 2014 2015 2016 
May [34.4,840.6] [47.875,3115.6] [39.06,3114.7] [9.56,287.6] [28.5,3115.6] 
June [ 81.7,865.5] [37.6,1606.5] [87.625,3115.6] [10.7,282.7] [30.687,394.2] 
July [197.3,6240.6] [40.375, 1903.1] [3111.9,3115.6] [6.56,3115.6] [35.125,1555.9] 

August [62.5,909.8] [175.1,2220.9] [81.375,3115.6] [550.56,3025] [57.625,3114.7] 
FCik [193-9,2214.1] [75.2,2211.5] [829.98,2336.5] (144.345, 1677.7] [37.98,2045.1] 


Particulate reduction-with an efficiency of 75%-and relative reduction ef- 
fect impact 0.2 and efficiency 40%, and the relative reduction effect impact 
0.375 have been shown by a green belt. 04, = [45,55] denotes environmen- 
tal loading capacities at area p in period k; transfer factor tp from source 
to area p is 1.1587 x 1072. The minimum and maximum wind speeds are 
[6 m/s, 15 m/s] during the 120 days; e, = [55,154] shows the PM10 emis- 
sion allowance during the period & (t/day) in Sistan. 


Table 10: PM10 efficiencies of different control measures 


Methods efficiency Interval efficiency 
Mulching method 7 (0.4, 0.6] 
Green belt method 12 (0.40, 0.75] 


We now solve the air pollution control model using the methods mentioned 
in Section 2. The results are given in Table 11. 


Table 11: Results obtained from solving ILP model (4) 


BWC ITSM UTSM 
ZX 10" | [0.1490,6.0239] | [2.19050165304,5.16560105314] | [2.50191551199,7.7566276346] 
Biopt 0.0,3348.7 [3634.1,3834.1] [3834.1 , 6795.9] 
D2opt 0.0,3463.2 [3730.7, 4830.7] [4830.7 , 8152.3] 
D3opt [136.0,4963.8] [4988, 5348] [5071.4 , 8480] 
Laopt 0.0, 2540.1 [3307.5,4147.5] [4147.5 , 7222.4] 
Deopt 0.0,3234.3 [3317.8,4117.8] [4117.8 , 7181.9] 
Zeopt 150.3,194.0 [193.95,663.7] [240.6 , 240.6] 
Dropt 75.2,120.4 [75.2,521.3] [46.8 , 46.8] 
Dsopt 20.9,1192.0 [16.2,447.4] [1.1044e-11 , 2.7188¢-07] 
Toopt 144,3,231.0 [144.4,618.9] [179.6 , 179.6] 
Fioopt 379.0,60.8 [3798,623.2] [185.4 , 185.4] 
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Results and analysis 

Tables related to the reduction of particulates are given in the previous sec- 
tion, using two controllers as mulching and green belt. It can be said that the 
Sistan region aerosols were considered over four months in five-year periods, 
such that the lower bound is the lowest amount of the particulate matter 
PM10 and upper bound represent the most amount of the particulate matter 
PM10 in Sistan. Now, using the two above mentioned controller, the reduc- 
tion of aerosols is calculated according to efficiencies in every method and 
these estimations are presented in Tables 8-11. Notably, the costs associated 
with each method can be considered based on the funding in the environ- 
mental section and different covering lands. Inflation in the country has been 
considered over the years according to variable data of costs and inflations 
as well as other parameters studied in the Sistan region. The greatest re- 
ductions in the particulate matter were calculated by the ILP model and the 
presented formulas. 


For greater reliability and accuracy of calculations, the ILP would be 
solved by three methods: BWC, ITSM, and UTSM. The variations in the 
reduction of aerosols are studied. The results show that if the aim is to min- 
imize the cost for decision-makers and officials, then the BWC method could 
be the best method because it can calculate the minimum and maximum 
of the costs under the two optimistic and pessimistic models (the best and 
worst) for finding the best answer. For example, the cost per hectare—using 
mulching-the minimum cost is estimated ~ 0.1490 x 107 and the maximum 
cost is estimated ~ 6.0239 x 107; Therefore, these are provided as an interval 
under the upper and lower bounds. 


In the BWC method, as mentioned in Section 2, the resulting solution 
does not exactly hold into the constraints. For example, in the first con- 
straint, the amount of reduction in particulates in every period may be more 
than the greatest reduction in particulates in each period. Therefore, UTSM 
and ITSM methods are used. In this case, maybe because of the reduction 
area which is created by the solutions, the model cannot be considered good 
enough in order to estimate the reduction of costs. It can make changes in 
the costs in some parts to find the best answer and increasing them because 
the constraints are not violated. Hence, we can say to find the best answer 
for the reduction of particulate matter, ITSM and UTSM methods are more 
suitable. In Table 11, 10 optimal solutions are given which are shown with the 
symbol x,,. The first five numbers represent 2,; according to the mulching 
and the second set of five numbers is 2 p,¢ according to the green belt over 
five years. As seen in the ITSM method, the changes are: [3634.1, 3834.1] in 
2012; in that year, the minimum value is 55 and the maximum value is 9985, 
which show the most amount of particulate matter in July. The above models 
considered all solutions in a feasible area for one another and the maximum 
value is 9985, which estimated the changes in the reduction of particulates 
(3634.14, 3834.14]. 
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ITSM method enabled shows a reduction in the intervals. Using the 
UTSM method, which can be said to be one of the best methods, it seeks 
to find the best solution in the feasible area. Optimally, the reduction of 
particulates is [3834.14, 6795.9128], which has obtained good results, and 
further reductions are anticipated. Also, in 2013, the minimum value is 62.5 
for the entire year and the maximum value is 4985, such that the reduction of 
the particulate matter has been estimated at [3730.7, 4830.7]. This indicates 
that in 2013, according to the amount of particulates, [3730.7, 4830.7] can 
undergo a decrease. It can be said that the obtained values from the model 
show a strong reaction over the years, according to the highest amount of 
the particulate matter and using the above analysis, we have the amount of 
particulates [4830.7 , 8152.3416] in the UTSM method similarly. It seeks to 
solve and attempts to reduce it, but according to the obtained numbers, the 
reduction is less in some years, it is because of the amount of particulate 
matter which is very large in all months of the year. However, there are the 
powerful and useful methods to control particulate matter, those may not 
be able to greatly reduce particulates, and many environmental factors can 
cause fewer reductions. The other numbers are explained similarly. 


4 Conclusion 


In this paper, some methods for solving the ILP models have been reviewed, 
such as the BWC, TSM, and improved TSM (ITSM). Besides, an updated 
two-step method (UTSM) has been proposed. This method improves the 
TSM by considering two feasibility and optimality conditions. The solutions 
comparison for the BWC, TSM, ITSM, and UTSM methods based on the 
feasibility and optimality conditions shows that some points of the BWC and 
TSM may be infeasible. Also, some solutions of the ITSM method may be 
non-optimal. The advantage of the UTSM is that, two constraints have been 
added in the second submodel which guarantees the feasibility and optimality 
conditions, and so the obtained region by the UTSM will be feasible and 
optimal. Also, a case study related to air quality management has been 
done. 

For solving many uncertain problems in the real-world, different models are 
used. The ILP model which is one of the models under uncertainty can be 
used for modeling these problems. In this paper, we consider an ILP model 
related to minimizing the costs in order to control aerosols and solve it by 
the BWC, ITSM, and UTSM methods. The results show that the lower and 
upper bounds of the costs obtain by the BWC method. The reduction of 
particulates can be clearly seen in two ITSM and UTSM methods. 
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